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We discuss multipole properties of filled skutterudites containing heavy lanthanide Ln from a micro- 
scopic viewpoint on the basis of a seven-orbital Anderson model. For Ln=Gd, in contrast to naive ex- 
pectation, quadrupole moments remain in addition to main dipole ones. For Ln=Ho, we find an exotic 
state governed by octupole moment. For Ln=Tb and Tm, no significant multipole moments appear at low 
temperatures, while for Ln=Dy, Er, and Yb, dipole and higher-order multipoles are dominant. We briefly 
discuss possible relevance of these multipole states with actual materials. 
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Recently, filled skutterudite compounds LnT4Xi2 with lan- 
thanide Ln, transition metal atom T, and pnictogen X have 
attracted much attention in the research field of condensed 
matter physics due to large varieties of exotic phenomena 
concerning magnetism and superconductivity. 1,2 Since these 
compounds crystallize in the cubic structure with high sym- 
metry of Th point group, 3 /-orbital degeneracy remains in 
general. Due to the strong spin-orbit coupling of / electrons, 
spin-orbital complex degrees of freedom, i.e., multipoles, be- 
come active in filled skutterudite compounds, which is a char- 
acteristic issue of this material group. 

In fact, recent experiments in cooperation with phenomeno- 
logical theory have revealed that multipoles appear ubiqui- 
tously in filled skutterudites. For instance, a rich phase di- 
agram of PrOs4Sbi2 including field-induced quadrupole or- 
der has been unveiled. 4-6 In NdFe4Pi2, a significant role of 
quadrupole at low temperatures has been suggested from the 
measurement of elastic constant. 7 In SmRu4Pi2, a possibility 
of octupole order has been proposed from experimental re- 
sults. 8-11 Recently, antiferro Ti-type order 12 has been actively 
discussed experimentally and theoretically for P1-RU4P12 13 ' 14 
and PrFe 4 Pi 2 . 15-17 

In such current circumstances, it is one of challenging is- 
sues to clarify microscopic aspects of multipole degrees of 
freedom in /-electron systems. The present author and col- 
laborators have developed a microscopic theory for multipole- 
related phenomena on the basis of a j-j coupling scheme. 18 In 
particular, possible multipole states for light lanthanide filled 
skutterudites have been discussed by analyzing multipole sus- 
ceptibility in detail. 19-25 

However, multipole phenomena in heavy lanthanide sys- 
tems have not been understood satisfactorily, mainly due to 
the difficulty in the theoretical treatment of multipoles in lan- 
thanides with large total angular momentum. Experimentally, 
due to the development of crystal growth technique under 
high pressure, heavy lanthanide filled skutterudite compounds 
have been synthesized 26-28 and research activities to clarify 
multipoles in such materials will be advanced. It is required 
to make more effort in the theoretical side to develop a micro- 
scopic theory for multipole state of heavy lanthanide systems. 

In this Letter, first we define multipole operators of / elec- 
trons. Then, we analyze multipole susceptibility of a seven- 
orbital Anderson model by using a numerical renormalization 
group technique. It is found that quadrupole moments appear 



even for Ln=Gd and an exotic phase governed by octupole 
moment is suggested for Ln=Ho. We find no significant mul- 
tipole moments at low temperatures for Ln=Tb and Tm, while 
dipole and higher-order multipole moments are dominant for 
Ln=Dy, Er, and Yb. We discuss these multipole states in com- 
parison with available experimental results. 
The seven-orbital Anderson model is given by 

H = £ kci a c ka + Y (V m c{j m « + h-c.) 
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where e k denotes conduction electron dispersion, c ka 
indicates the annihilation operator for conduction elec- 
tron with momentum k and spin a, a=+l (—1) for up 
(down) spin, f ma is the annihilation operator for / elec- 
tron with spin a and z-component m of angular mo- 
mentum £=3, V m is the hybridization between conduc- 
tion and / electrons, £> m . m < is a crystalline electric field 
(CEF) potential for / electrons, 8 aa i is the Kronecker's 
delta, A is the spin-orbit coupling, Cm,±i,m,±i=±w/2, 
Cm±i,=Fi,m,±i=\/l2 — m(m ± l)/2, and zero for the other 
cases. The Coulomb integral I is expressed by the combina- 
tion of Slater-Condon parameters, F°, F 2 , F 4 , and F 6 . 29 

For filled skutterudites, the main conduction band is given 
by a u , constructed from p-orbitals of pnictogen. 30 Note that 
the hybridization occurs between the states with the same 
symmetry. Since the a u conduction band has xyz symme- 
try, we set V2=—V-2=V and zero for other to. For the Th 
point group, B mm / is given by three CEF parameters, B®, 
£g,and B 2 , 3 expressed as B%=Wx/15, Bg=W(l-|x|)/180, 
and S^W^y/24. 31 ' 32 Here x and y specify the CEF scheme, 
while W determines the energy scale of the CEF potential. In 
order to adjust the local /-electron number n, we appropri- 
ately change the chemical potential. 

Let us set the parameters of the model. Details of the deter- 
mination of the parameters will be discussed elsewhere. Con- 
cerning Slater-Condon parameters, we set F°=10 eV by hand. 
Others are determined so as to reproduce excitation spectra 



.3+ 33 

F 4 =6.60 eV, and F 6 =4.44 eV. We use these values for all 



of Pr J+ ion. Here we show only the results: F =8.75 eV, 
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lanthanides. On the other hand, we use experimental value 
of A for each lanthanide such as A=0.180 eV (Gd), 0.212 eV 
(Tb), 0.240 eV (Dy), 0.265 eV (Ho), 0.295 eV (Er), 0.326 
eV (Tm), and 0.356 eV (Yb). 34 The hybridization is fixed as 
y=0.05 eV and a half of the bandwidth of a u conduction band 
is set as 1 eV. Concerning CEF parameters, we set W=— 0.4 
meV, y=0.3, and x=0.3 so as to reproduce quasi-quartet CEF 
scheme of PrOs 4 Sb 12 . 3 " 7 

Here we provide a comment on the local CEF state. Due 
to the lack of space, we do not show the CEF energy schemes 
for n=7<~13, but they are almost the same as those obtained in 
Ref. 20, even though we have used different parameters there. 
Namely, the CEF energy schemes are not changed drastically, 
as long as we use realistic Coulomb interaction and spin-orbit 
coupling. As for the CEF ground state, except for n=10, it 
is not easily converted when we slightly change the value of 
x around at £=0.3. Thus, we set £=0.3 as a typical value for 
n=7^9 and 1 1~13. For n=10, it has been found that r^ 3 dou- 
blet and r+ triplet states are almost degenerate in the wide 
range of the values of x. 20 Namely, it is necessary to discuss 
carefully the case of n=10. 

Now we define the multipole operator X for / electron. In 
general, X is expressed as 

X = ^>«Tf\ (2) 

where k is a rank of multipole, 7 is a label to express Oh 
irreducible representation, and is cubic tensor operator, 
given by Tj = \, q G\, q T q . Here an integer q runs between 

" (k) (k) 

—k and k, T q is spherical tensor operator, and G\ t ' q is the 
transformation matrix between spherical and cubic harmon- 
ics. The coefficient p^ will be discussed later. 

In this paper, we define multipole as spin-orbital density in 
the form of one-body operator from the viewpoint of multi- 
pole expansion of electron density in electromagnetism. Then, 

" (k) 

the spherical tensor operator T q for / electron is given in the 
second-quantized form as 

f(k) _ v T {k ' q) ft f / , (3) 

± q ~ / j - L mtT,m'i7' JmaJm'cr' , WJ 

(k a) 

where T^' ma , can be calculated by using the Wigner-Eckart 
theorem as 38 

x {jn\ems-){jii'\£m's—). 

Here £=3, s=V2, j=£±s, fi is the z-component of j, 
(j lAj' j" V>") denotes the Clebsch-Gordan coeffi- 
cient, and is the reduced matrix element 
for spherical tensor operator, given by (j\ \ \j)= 
y/(2j + k + l)!/(2j - k)\/2 k . 3S Note that k<2j and the 
highest rank of /-electron multipole is 7 in the present 
definition. 39 

We note that it is not necessary to take double summations 
concerning j in eq. (4), since the matrix representation of to- 
tal angular momentum J is block-diagonalized in the (j, y)- 
basis. We have checked that the same results as eq. (4) are ob- 
tained when we calculate higher-order multipole operators by 



following the symmetrized expression of multiple products of 
^38,40,41 given in ^ ( TO; o-^basis as J TOCTiTOV /= 5„„>L mm >+ 
Saa' 5mm 1 ■ Here L denotes angular momentum operator for 
£=3 and S indicates spin operator, given by S=er/2 with Pauli 
matrix er. 

(k) 

Let us now determine the coefficient p\ ' . In order to dis- 
cuss the multipole state, it is necessary to evaluate the mul- 
tipole susceptibility in the linear response theory. 21 However, 
multipoles belonging to the same symmetry are mixed in gen- 
eral, even if the rank is different. In addition, multipoles are 

(k) 

also mixed due to the CEF effect. Thus, we determine p\ by 
the normalized eigenstate of susceptibility matrix 

1 p-Ei/T _ p-Ej/T 

Xw4E f F (i\m k) ~ P { 7 k) }\3) 

A > Ej ~ Ei (5) 

where Ei is the eigenenergy for the z-th eigenstate \i) of H, 

T is a temperature, p {k) =Y,i e~ E ' /T (i\f( k) \i)/Z, and Z is 
the partition function given by Z=]T\ e - £ */ T . Note that the 
multipole susceptibility is given by the eigenvalue of the sus- 
ceptibility matrix. 

In order to evaluate the susceptibility matrix, here we em- 
ploy a numerical renormalization group (NRG) method, 42 
in which momentum space is logarithmically discretized to 
include efficiently the conduction electrons near the Fermi 
energy. In actual calculations, we introduce a cut-off A for 
the logarithmic discretization of the conduction band. Due to 
the limitation of computer resources, we keep M low-energy 
states. In this paper, we set A=5 and M=4500. Note that the 
temperature T is defined as r=A-( Ar - 1 )/ 2 in the NRG calcu- 
lation, where N is the number of the renormalization step. 

Here five comments are in order, (i) In this paper, we 
call higher multipoles by following the numerical terms of 
IUPAC, 43 i.e., 16=hexadeca, 32=dotriaconta, 64=tetrahexa- 
conta, and 128=octacosahecta. (ii) We express the irreducible 
representation of the CEF state by Bethe notation in this pa- 
per, but for multipoles, we use short-hand notations by the 
combination of the number of irreducible representation and 
the parity of time reversal symmetry, g for gerade and u for 
ungerade. (iii) For Th symmetry, T\ and T 2 of Oh are mixed. 
We remark that T4 and T5 of Oh are also mixed. 3 Thus, we 
obtain six independent multipole components as lg+2g, 2u, 
3g, 3u, 4g+5g, and 4u+5u for filled skutterudites. Note that 
lu does not appear within rank 7. (iv) We normalize each 
multipole operator so as to satisfy the orthonormal condi- 
tion Tr{T^Ty ''}= dkk'S-yy, 44 when we express the mul- 
tipole moment as eq. (2). (v) The susceptibility for 4u mul- 
tipole moment does not mean magnetic susceptibility, which 
is evaluated by the response of magnetic moment L+2S, i.e., 
J+S. 19 ' 20 

Now we discuss the numerical results. In Fig. 1 (a), we show 
the results for n=7 (Gd 3+ ion). Note that we depict only the 
largest eigenvalue state of each symmetry group. The CEF 
ground state is well described by L=0 and J=S=1I2, but it 
is almost independent of x, since the CEF potentials for L=0 
provide only the energy shift. The dominant multipole com- 
ponent is 4u dipole and the contributions from higher multi- 
poles are negligibly small. We find that the secondary com- 
ponents are 3g and 5g quadrupoles. Around at T=10 -7 , there 
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Fig. 1. (Color online) T\ for (a) n=7, (b) n=8, and (c) n=9. 



occurs partial screening due to conduction electrons and at 
low temperatures, spin S=l/2 remains. 19,20 In such a region, 
quadrupole components disappear. 

In the LS coupling scheme, the / 7 state is uniquely speci- 
fied by J=S=l/2 and L=0. Thus, quadrupole does not seem to 
appear at the first glance. However, we should note that actual 
situation is always deviated from the LS coupling scheme. 23 
Some finite contribution of the j-j coupling scheme is in- 
cluded in the ground state. If we consider very large A, first 
we fully occupy j=5/2 sextet and then, we accommodate one 
/ electron in j=l/2 octet. Thus, this state can be multipole- 
active. Quite recently, it has been observed in GdRu4P 12 that 
101 Ru NQR frequency exhibits temperature dependence be- 
low a Neel temperature Tn=22K. 45 This may be interpreted 
as the effect of quadrupole due to the deviation from the LS 
coupling scheme. 

In Fig. 1(b), we show the results for n=8 (Tb 3+ ion). At 
x=0.3, the CEF ground state is Tf singlet and the excited state 
is triplet with excitation energy of 1 meV. For T>10~ 4 , 
the dominant component is 4u dipole, while for T<10~ 4 , we 




1 °1 oVioVio^ioWio^io 



Fig. 2. (Color online) T\ for n=10 with (a) x=0.3 and (b) x=0.25. 



find no significant local moments. Note that we obtain finite 
values for p^ , pf^ , and pf^ due to the effect of CEF poten- 
tials. Thus, we cannot exclude a possibility of antiferro lg+2g 
ordering. 1316 For TbRu 4 Pi 2 , T N is found to be 20K, but an- 
other anomaly is observed at Ti = 10K. 26 It may be interest- 
ing to discuss Ti in the context of antiferro hexadecapole and 
tetrahexacontapole order. 

In Fig. 1(c), the results for n=9 (Dy 3+ ion) are depicted. In 
this case, the CEF ground state is doublet and the excited 
state is quartet with excitation energy of 2.7 meV. The 
dominant multipole is always the mixture of 4u and 5u. The 
main component is dipole, but we find significant contribution 
from octacosahectapole. 

Next we focus on the case of n=10 (Ho 3+ ion). At x=0.3, 
the CEF ground state is triplet, but the first excited state 
is r^g doublet with very small excitation energy such as 10~ 5 
eV. Thus, Ho-based filled skutterudite is considered to be in 
the quasi-quintet situation. 20 ' 46 In Fig. 2(a), we show the re- 
sults for n=10 and £=0.3. We observe several kinds of mul- 
tipoles, but the dominant one is always given by the mixture 
of 4u and 5u from dipole, octupole, dotriacontapole, and oc- 
tacosahectapole. 

As mentioned above, the CEF ground state is fragile for 
n=10. If we slightly decrease x, the ground state is eas- 
ily changed. Then, we evaluate multipole susceptibility for 
£=0.25 with doublet ground state, as shown in Fig. 2(b). 
The multipole states at high temperatures are similar to those 
in Fig. 2(a), but at low temperatures, we find 2u multipole 

(3) (7) 

state, expressed as p 2u =0.955 and p 2u =0.297. Namely, the 
main component is 2u octupole, but there is contribution of 2u 
octacosahectapole. From the elastic constant measurement for 
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Fig. 3. (Color online) T\ for (a) n=l 1, (b) n=12, and (c) n=13. 



HoFe4Pi2, some anomalous features have been discussed. 46 It 
may be interesting to reexamine experimental results from the 
viewpoint of 2u octupole state. 

Let us move onto the cases of n=l 1, 12, and 13. In Fig. 3(a), 
we show the results for n=l 1 (Er 3+ ion). For n=\ 1, the CEF 
ground state is T§ doublet and the excited state is Tq 7 quartet. 
The excitation energy is 5.6 meV. We find that the dominant 
multipole state is given by 4u+5u, but main components are 
4u dipole and 4u octacosahectapole. The admixture is almost 
independent of temperature. 

In Fig. 3(b), the results for n=12 (Tm 3+ ion) are depicted. 
The CEF ground state at x=0.3 is singlet, while the ex- 
cited state is triplet with the excitation energy of 5.9 meV. 
In the high-temperature region, several kinds of multipoles 
appear, but when we decrease temperature, there appear no 
significant multipole moments, as in the case of n=8. Again 
we mention that a possibility of antiferro lg+2g ordering can- 
not be excluded. 

In Fig. 3(c), we depict the results for ?i=13 (Yb 3+ ion). 
The CEF ground and excited states are both T§ doublets with 



the excitation energy of 5.1 meV. When temperature is de- 
creased, there appears the dominant mixed multipole of 4u 
and 5u with significant components of dipole, octupole, dotri- 
acontapole, and octacosahectapole. Among them, large con- 
tributions come from dotriacontapole and octacosahectapole. 

Finally, let us briefly discuss the Kondo effect. First we 
note that the so-called underscreening Kondo effect can oc- 
cur in the present model, since we include only the single a u 
conduction band hybridized with some of seven / orbitals. 
In fact, we have found the partial screening for rv=l, but for 
other values of n, we did not observe the Kondo behavior in 
the present temperature range, since the Kondo temperature 
becomes much smaller than the CEF splitting. However, if V 
is increased, we expect to find the Kondo behavior for n^7. 
When we perform the calculations for V=0.5, we actually ob- 
serve that the residual 4u moment is eventually screened even 
in the present temperature range for n^=l . Note, however, that 
the Kondo behavior can be found only for unrealistic large 
value of V in the present model. In order to discuss appro- 
priately the Kondo effect in filled skutterudites, it is highly 
required to include not only a u but also e u conduction bands. 
It is one of future problems. 

In summary, we have discussed possible multipole states of 
filled skutterudites including heavy lanthanide ion Ln 3+ . For 
Ln=Gd, we have found the effect of quadrupole moments due 
to the deviation from the LS coupling scheme. For Ln=Ho, 
the CEF ground state is quasi-quintet with doublet and 
triplet. When we slightly change the CEF potential, the exotic 
state dominated by 2u octupole moment has been found. For 
Ln=Tb and Tm, the CEF ground state is singlet and we 
have found no significant multipole moments at low temper- 
atures, although we cannot exclude a possibility of antiferro 
lg+2g ordering. For Ln=Dy, Er, and Yb, the CEF ground state 
is doublet, and the dominant moment is the mixture of 
4u and 5u. We believe that the present results stimulate fur- 
ther study of multipole phenomena in heavy lanthanide filled 
skutterudites. 
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by the Ministry of Education, Culture, Sports, Science, and 
Technology of Japan. The computation in this work has been 
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